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Abstract 

We present a Monte-Carlo simulation of energy deposition process in relativistic 
heavy-ion collisions based on a new realization of the Interacting-Gluon-Model (IGM) 
for high energy N — N collisions. In particular we show results for proton spectra 
from collisions of Eiab = 200 GeV/N '^^S beam incident on '^^S target and analyze the 
(^f^ energy and mass dependence of nuclear stopping power predicted by our model. Theo- 

retical predictions for proton rapidity distributions of both ^"^Pb -|- ^"^Pb collisions at 
Eiab = 160 GeV/N CERN SPS and ^^^Au + ^^^Au at = 200 GeV BNL RHIC 

are given. 

in 

^ PACS: 25.75.-Fr 

^ 1 Introduction 

Whether the central rapidity region is baryon free or not after the collision is an impor- 
tant issue in the field of relativistic heavy-ion collisions. These two different pictures imply 
different mechanisms for QGP formation, if QGP is really produced there, and thus will 
make the detection of QGP different for the two situations. The subject is directly related 
to how the C.M.S energies of the colliding nuclei are deposited in this region, i.e., to the 
nuclear stopping power. Though much work that has been done on both theoretical and 
experimental aspects of this subject indicates so far a baryon rich picture (i.e., complete 
stopping) [1] for heavy systems such as gold on gold below the AGS energies, where initial 
conditions for QGP formation are not well fulfilled, a common consensus on the issue has 
not yet been reached for heavy-ion collisions at CERN SPS to BNL RHIC and CERN LHC 
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energies, where one expects the formation of QGP. For example, FRITIOF [2] predicts fully 
baryon free in the central rapidity area for collisions of gold on gold at RHIC energies while 
VENUS [3] and DTU [4]) predict the opposite, although except for this, all models have 
almost the same predicting abilities from AGS to SPS energies. 

Whereas in the above mentioned models the main degrees of freedom are quarks (va- 
lence in LUND and with the sea-quark admixture in VENUS and DTU) and gluons are 
treated collectively as non-interacting strings, the Interacting-Gluon-Model (IGM) [5, 6, 7] 
pays more attention to the role of gluons in high energy collisions, following the fact that 
in the perturbative QCD (PQCD) gluon-gluon interactions are much more important than 
those between quarks [8]. IGM was already applied successfully to nuclear collisions, how- 
ever either in the coherent nuclear tube version which does not allow for calculation of the 
nucleonic spectra [6] or in the version devoted to the specific problem of looking for coherent 
vs chaotic aspects of nuclear collisions, also not directly suitable for our present purposes [7]. 

As both Ref.[6] and Ref.[7] have already shown the basic ability of IGM to describe 
adequately essential details of multiparticle production processes, we shall concentrate in 
this letter only on the study of nuclear stopping power provided by IGM when extended to 
heavy-ion collisions at C.M.S energies no smaller than those of SPS. 

2 Brief review on the IGM for high energy N — N pro- 
cess 

The IGM can be summarized as follows [5]: According to PQCD: agg ^ agq ^ aqq, where 
agg, agq aud aqq are the interaction cross sections of a gluon-gluon, gluon-quark and quark- 
quark pair, respectively. Therefore, assuming that the above inequality holds also in soft 
interactions, in each event one has he following picture [5]: 

• Two colliding hadrons are represented by the valence quarks carrying their quantum 
numbers plus the accompanying clouds of gluons (which also include the sea qq pairs 
and should therefore be regarded as effective ones). In the course of the collision, 
gluonic clouds interact strongly and form rmm-firehalls (MF's) which are supposed to 
populate mainly in the central rapidity region of the reaction. 

• The valence quarks (together with the gluons which did not interact) get excited and 
form leading jets (LJ's) (or beam jets) which are supposed to populate mainly in the 
fragmentation regions of the reaction. Here we shall assume further that LJ's are 
identical with leading particles (practically nucleons). 

The above mentioned MF's eventually form a lump of gluonic matter which is called a cen- 
tral fireball (CF). It is this CF that governs the final state particle production in the central 
rapidity region. 

In the framework of the IGM, if fractions of x and y of the energy momentum of the 
incoming nucleons are deposited in central rapidity region and form a CF, then the energy 
E and momentum P and invariant mass M and rapidity A of this CF are expressed as (all 
masses are neglected here): 




(1) 



and 





(2) 
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while the probability distribution for producing such a CF is [6, 11] 



X{x, y) = xo ^M--^^ - ] dixy - KLu) (3) 

where xo is a normalization factor and Kmin = Mq / ^ stands for the minimal inelasticity 
allowed. There is simple relation among the moments [x) , (x^) and (y), (j/^), which is valid 
at high energy, namely 

(.) c 2{x') c (y) c 2{y') = J^^^^^ (4) 

where a is a constant related to the gluonic cross-section, ph represents fraction of the 
hadronic energy-momentum attributed to gluons in hadron structure function [9] (in fact 
we have here p'f^ = PprojectUe • Ptarget = p^ ' P^), ^0 is the lightest mass of the possible CF 
and (rf^^ (s) is the inelastic cross section for high energy N — N process. 

We have thus essentially three input parameters (although highly correlated, cf. Eq.(4)): 
Mq delineating the gluonic phase space (i.e., x y > K^i„ = Mq/s), a fixing the strength of 
gluonic interaction and product of p^'s dictating what part of energy-momentum of initial 
nucleons is used to form CF. ^ 



3 IGM and high energy heavy-ion collisions 

We proceed now to simulation of energy deposition process in relativistic collisions of two 
nuclei A and B (with nucleus B incident on nucleus A), where A and B denote also their 
corresponding atomic numbers. In IGM it is realized in steps according to the following 
algorithm [13]: 

(i) We simulate positions of all the A and B nucleons, in the rest frames of colliding nuclei 
with the respective coordinate origins located at their mass centers and according to their 
nucleon densities Pa,b(^) given by the usual three-parameter Woods-Saxon form [14]. 

(ii) For every heavy-ion collision event we sample the impact parameter b using the 
Glauber model [15], in which the inelastic cross section of high energy A — B process is 
given by 

dh[l- f-'^^n^^ABih)^ (5) 



with 



Tab{'^) = j dhji^ dzA dh-Q dzB PAit>ji^,ZA) PBit>-Q,ZB) 6{h -hji^ + h-Q) (6) 

being the thickness function normalized to 1. Consequently, the probability distribution for 
the impact parameter b is given by 



dP{h) 1-e-^in TABih) 

- -AR • (7) 



dh 



^It is perhaps worth mentioning that such a formulated IGM was, in addition to Refs.[[5, 6, 7]], used 
also to discuss the problem of inelasticity in cosmic ray experiments [10], to connect the limited number of 
clans in multiparticle distributions with the limited behaviour of inelasticity [11], and recently it was also 
discussed in connection with the possible mini-jets production [12]. 
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Once the impact parameter b is selected according to Eq.(7), the transverse and longitu- 
dinal coordinates of all colliding nucleons can be determined from their initial coordinates. 
Finally all nucleons in both colliding nuclei are ordered according to their z-coordinates in 
lab. frame, which results in: zp(B) < zp(B — 1) < zp(B — 2) < • • • < zp(l) < zt(1) < 
ztC^) < ■ ■ ■ < zt(A) (zp(i) denotes the z-coordinate of the ith nucleon in nucleus B with 
i G (1, B) and zpij) stands for the z-coordinate of jth nucleon in nucleus A with j G (1, A)). 

(iii) Finally, we simulate the nucleus-nucleus collision as proceeding via a number of 
consecutive N — N like collisions in the following way: 

• For a pair of nucleons ij: i = 1, . . . , B from nucleus B and j = 1, . . .A from nucleus 
A, we check first if they interact. They do it if their transverse distance is smaller 

I ^NN (g\ 

than Y — - and -^5^" > 2.14 GeV . Otherwise this pair will not collide and the 
four-momenta of each of the two nucleons are kept unchanged. Here ^Js^ stands for the 
energy of the k-ih pair in its cm. frame. Both impinging nucleons and the produced 
leading particles are assumed to propagate along straight lines. 

• For the pair of nucleons which do interact, we carry out the corresponding nucleon- 
nucleon collision which proceed at current energy of the colliding nucleonic pair, 

whether th collision is elastic or inelastic is sampled according to the ratio 
of the corresponding cross sections a^^^ {sj.) and a^J^ {sj.) taken from experimen- 
tal data [16]. For inelastic collision, the IGM is called for providing a new gluonic 
CF and two new leading particles formed from both participating nucleons. The 
energy-momentum conservation implies the continuous degradation of the collision 
energy of interacting nucleonic pairs along the chain of consecutive scatterings. As 
all the 3 essential parameters of the IGM are kept unchanged with k (and equal to: 
ap\ = 0.05 GeV^ fm'^ , Mq = 0.22 GeV)^ , the above procedure corresponds to a tacit 
assumption of restoration of the gluonic clouds in the nucleons propagating through 
the nuclear matter after every inelastic collision. The magnitudes of the corresponding 
transverse momenta of colliding nucleons are sampled according to a Gaussian distri- 
bution with (pt) = 0.45 GeV/c [17]. For elastic scattering, we use differential cross 
section from Ref.[18] to generate new four-momenta of colliding nucleons. 

Such an implementation of IGM to nuclear collisions differs distinctly from that in Ref.[6] 
or Ref.[7]. In fact, it can be regarded as a new realization of the IGM, similar to the so 
called "model C"[19]. Contrary to Refs.[[6, 7]], it allows the presence of additional energy- 
momentum transfer from the valence quark component of the colliding nucleons to gluons 
during the nuclear collision process (specifically: in between the consecutive scatterings). 
As we shall see below, such a scenario leads to the nuclear stopping compatible with exper- 
imental data. 



4 Results and discussion 

Fig.l and Fig. 2 display the event-normalized proton rapidity distributions in central colli- 
sions (corresponding to zero impact parameter, b = 0, in the above IGM.) of ^^S -|- ^^S and 
13"^ Au + 13"^ Au, at CERN SPS and BNL RHIC energies, respectively. The noticeable feature 
of the displayed results is the increase of baryon number in the central rapidity region with 
increasing atomic mass of colliding nuclei and the decrease of baryon number in the region 
with growing energy. 



^As mentioned before, a and p/^'s are strongly correlated and therefore usually given as product like in 
this case above. If we assume that like in Ref.[5], p/j ~ 0.5, then a = 0.2 GeV^ fm^ . 
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A useful measure [20, 21] of the nuclear stoping power characterizing the baryon content 
of the central rapidity region is the mean rapidity shift (Aj/p) of the projectile participant 
protons from their original beam rapidity. The larger is this value, the stronger is the 
stopping of the colliding system and the more baryons are showing up in the central rapidity 
region. In order to compare this measurement to those of other colliding systems at different 
energies, a relative stopping parameter, which also is an alternative of the event-normalized 
proton rapidity distributions, was introduced as follows[22]: 



Vbeam Vcm 

with j/cm being the rapidity of center of mass of all participant nucleons. The values of S for 
collisions of two mass symmetric nuclei as calculated in our model at CERN SPS and BNL 
RHIC energies are tabulated in Table 1. They confirm the statements suggested already by 
Fig.l and Fig. 2, namely that heavy systems are more stopped than light ones at the same 
energy -yiiviv whereas for the same colliding system, the larger the -yiiviv, the weaker is the 
nuclear stopping power. 



Table 1: The energy and mass dependence of relative stopping parameter S. 



System 




b(fm) 


Theor. 


Exp. 


32S + 32S 


20 


0-1 


0.51 


0.53 


32S + 3^S 


20 





0.49 


No 




20 





0.67 


No 


32S + 32S 


200 





0.35 


No 


197^^^197^^ 


200 





0.59 


No 



The first conjecture results from the successive multiple scatterings of colliding nucleons 
and this effect manifested itself in recent experimental data [21] from NA35 Collaboration. 
The effect would increase the number of collisions an incident nucleon experiences, and 
would become stronger for heavier colliding systems. Therefore with a prion more inelastic 
collisions (which means more energy loss) present, heavy system are stopped more than 
light one at the same colliding energy. In fixed target Eiab = 200 GeV/N ^^S + ^^S central 
collisions, there is good agreement between theoretical calculations and experimental results 
both for the relative stopping parameter S as shown in Table. 1 and for the proton spectra 
as shown in Fig. 3 and Fig. 4. This agreement indicates that the multiple scattering effect 
has been properly accounted for in our simulation. 

The second statement is caused by the energy dependence of cross section (rf^^ (s) in 
elementary high energy N — N processes. The smaller energy loss of a participant projectile 
proton or a weaker stopping power corresponds to a smaller value of (x), and vice versa. 
Therefore a larger energy y^s^N of the colliding system leading to a larger (rf^^ (s) according 

2 

to Ref.[16], would results in a smaller value of (x) according to Eq.(4): (x) ~ jf2'^NNjj'y and 
thus leads to a weaker nuclear stopping power. 

As a theoretical prediction we shown in Fig. 5 the event-normalized proton rapidity dis- 
tributions in central collisions of Eiah = 160 GeV/N ^°^Fh + ^°^Fh with hfh = fm, for 
which we predict the relative stopping parameter S = 0.67. 
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The theoretical results concerning proton spectra shown in Fig. 3 and Fig. 4 are both 
normalized to the experimentally observed proton number in order to facilitate comparison 
with data. Because in this letter we have assumed for simplicity that all participant protons 
are observed in the final state of the heavy-ion collisions, our estimations provide only upper 
limits for proton rapidity distributions. However, in mass symmetric heavy-ion collisions, it 
is expected to result in a reasonable conclusion on whether or not the central rapidity re- 
gion is baryon free. This happens because other baryon productions (mainly strange baryon 
productions) which contribute to the net baryon rapidity distributions, proceed via fragmen- 
tation of the excited forward (backward) baryonic strings, (cf., Ref.[3]), and thus manifested 
themselves mostly in the projectile (target) rapidity region (also their abundances are much 
less than that of the protons considered there). ^ 



5 Conclusions 

Extending the IGM for high energy N — N process to relativistic heavy-ion collisions, we 
simulated energy deposition process in high energy nucleus-nucleus collisions. Our theoret- 
ical results of proton spectra in Eiab = 200 GeV/N ^^S -|- ^^S central collisions are in good 
agreement with experimental data. Moreover, according to our calculations, the amount of 
baryons present in the central rapidity region in relativistic heavy-ion collisions, is deter- 
mined both by the mass and by the energy -yiiviv of the colliding systems. It decreases 
with the increase of -yiiviv (indicating a weaker nuclear stopping power) and dramatically 
increases with the mass of the colliding systems (the final state baryons are inclined to be 
piled up in the central rapidity region signifying a stronger nuclear stopping power). Finally, 
if QGP will really be produced in Eiah = 160 GeV/N ^°^Fh + ^°^Fh central collisions at 
SPS or in ^iiviv = 200 GeV ^^^Au -|- ^^^Au central collisions at RHIC, it is predicted 
to be a baryon rich one, not a baryon free one. Although in this respect our conclusion is 
generally similar to that of the string models (cf., VENUS [3], DTU [4] and also [23]), we 
foresee a difference between IGM and the above mentioned models in the LHC energy range. 
The reason is the limited (i.e., at most increasing to a limited value) energy behaviour of 
inelasticity parameter of hadronic collisions predicted by IGM [10, 12] whereas it increases 
rather strongly in models presented in Refs.[[3, 4, 23]] (cf. Ref.[10] for more details). It 
means that asymptotically we should have generally broader distribution of baryons than 
those models predict, with width asymptotically constant (versus decreasing one in their 
case). 

Acknowledgement: This work is partly supported by Academia Sinica and National 
Science Foundation of China. 
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Figure captions 



Fig. 1 Theoretical results of event-normalized proton rapidity distributions for fixed target 
central collisions (& = fm) of Eiah = 200 GeV/N ^^S + ^^S (solid line) and ^^^Au + 
is^Au (dotted line). 

Fig. 2 Theoretical results of event-normalized proton rapidity distributions for central col- 
lisions (6 = fm) of ^^S -|- ^^S (shown in solid line) and "'^^^Au -|- "'^^^Au (shown in 
dotted line) collisions at RHIC energy -yiiviv = 200 GeV . 

Fig. 3 Event-normalized proton rapidity distributions in Eiab = 200 GeV/N ^^S -|- ^^S 
fixed target central collisions. Solid line stands for theoretical result normalized to 
the experimentally observed number of protons. Experimental data are from reference 
[21]. 

Fig. 4 Event-normalized proton transverse momentum distributions in fixed target central 
collisions of ^^S -|- ^^S with Eiab = 200 GeV/N. Solid line stands for theoretical result 
normalized to the experimentally observed number of protons. Experimental data are 
from reference [21]. 

Fig. 5 Theoretical prediction for event-normalized proton rapidity distributions in Eiah = 
160 GeV/N and & = fm fixed target ^°^Fh + ^°^Fh collisions. 
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